Superconducting d-wave stripes in cuprates: 
Valence bond order coexisting with nodal quasiparticles 
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We point out that unidirectional bond-centered charge-density-wave states in cuprates involve 
electronic order in both s- and d-wave channels, with non-local Coulomb repulsion suppressing the 
s-wave component. The resulting bond-charge-density wave, coexisting with superconductivity, is 
compatible with recent photoemission and tunneling data and as well as neutron-scattering measure- 
ments, once long-range order is destroyed by slow fluctuations or glassy disorder. In particular, the 
real-space structure of d-wave stripes is consistent with the scanning-tunneling-microscopy measure- 
ments on both underdoped Bi2Sr2CaCu20g+i and Ca2-iNa I Cu02Cl2 of Kohsaka et al. [Science 
315, 1380 (2007)]. 
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I. INTRODUCTION 

A remarkable aspect of the copper-oxide high-T c su- 
perconductors is that various ordering phenomena appar- 
ently compete, including commensurate and incommen- 
surate magnetism, superconducting pairing, and charge- 
density-wave formation. (More exotic states have also 
been proposed, but not verified experimentally beyond 
doubt.) While commensurate magnetism and supercon- 
ductivity are common phases in essentially all families 
of cuprates, the role of other instabilities for the global 
features of the phase diagram is less clear. 

A particularly interesting role is taken by charge- 
density waves. Such states break the discrete lat- 
tice translation symmetry, with examples being stripe, 
checkerboard, and valence-bond order. In the compounds 
La2-sBa K Cu04 and La2- x Sr :E Cu04 (with Nd or Eu co- 
doping) evidence for stripe-like spin and charge mod- 
ulations with static long-range order were detected, 1-4 
in particular near l/8th doping. (This is supported 
e.g. by strong phonon anomalies seen in neutron scat- 
tering experiments. ) While in other cuprate families 
similar long-range order has not been found, signa- 
tures of short-range charge order, likely pinned by im- 
purities, have been observed in scanning tunneling mi- 
croscopy (STM) on underdoped Bi2Sr2CaCu20s+5 6_9 
and Ca2_2;Na x Cu02Cl2. 9 ' 10 The low-energy electronic 
structure in the presence of charge order turns out 
to be remarkable: In La 15 / 8 Ba 1 /gCu04, angle-resolved 
photoemission (ARPES) indicated a quasiparticlc (QP) 
gap with d-wave like form, i.e., charge order coexists 
with gapless (nodal) QP in (1,1) direction (while antin- 
odal QP near (0, n) are gapped). 11 STM data on both 
underdoped Bi 2 Sr 2 CaCu208+5 and Ca2- :c Na :c Cu02Cl2 
show QP interference arising from coherent low-energy 
states near the nodes, whereas electronic states at 
higher energy and wavevectors close to the antinode are 
dominated by the real-space modulation of the short- 
range charge order. 7,9,12 This dichotomy in momen- 
tum space is has also been found in ARPES exper- 
iments in L^-zSr^CuC^, 13 Bi2Sr2CaCu20s+<5, 14 and 
Ca2_ 2; Na 2 ;Cu02Cl2 15 where well-defined nodal and ill- 



defined antinodal QP are frequently observed. 

These results suggest that momentum-space differenti- 
ation and tendencies toward charge ordering are common 
to underdoped cuprates. 16-18 The concept of fluctuating 
stripes, i.e., almost charge-ordered states, has been dis- 
cussed early on }^< 1G < ls > 19 This concept, appropriate for 
compounds without static long-range order, assumes the 
existence of a nearby stripe-ordered state, with physical 
observables being influenced by the low-lying collective 
modes associated with a charge-ordering instability. Fol- 
lowing this idea, we have recently calculated 20 the spin 
excitation spectrum of slowly fluctuating (or disordered) 
stripes. We were able to show that fluctuating stripes 
give rise to an "hour-glass" magnetic spectrum, very sim- 
ilar to that observed in neutron scattering experiments 
both on La2_ x Baa;Cu04 21 and YBa 2 Cu 3 6+ 5. 22,23 

The focus of this paper is on the electronic structure 
of stripe states. We introduce the concept of "ci-wave 
stripes" : Here the modulation of charge densities has pri- 
marily a d-wave form factor, i.e., lives more on the bonds 
than on the sites of the square lattice, leaving nodal QP 
unaffected. We illustrate that a picture of such bond- 
centered charge order, coexisting with superconductivity 
(this state may be dubbed "valence-bond supersolid") 
is consistent with various features seen in both ARPES 
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FIG. 1: Schematic real-space structure of a stripe state 
with primarily d-wave character and a 4 x 1 unit cell, i.e. 
Q = (±7r/2, 0). Cu lattice sites are shown as circles, with their 
size representing the on-site hole densities. The line strengths 
indicate the amplitude of bond variables like kinetic and mag- 
netic energies. The modulation in the site charge densities is 
small, whereas the one in the bond densities is large and of 
d-wave type. 24 Note the similarity of the bond modulation 
with the STM data of Ref. 9. 
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and STM measurements. In particular, the real-space 
pattern of (Z-wave stripes, Fig. 1, is strikingly similar 
to the STM results of Ref. 9, obtained on undcrdoped 
Bi 2 Sr 2 CaCu208+(5 and Ca 2 - x Na x Cu02Cl2. 

Other types of (Z-wave particle-hole order have been 
discussed before, (a) The (Z-density wave (or staggercd- 
fiux) phase 25 was proposed as a candidate for the pseudo- 
gap. However, it has no measurable charge modulation 
and shall not be discussed here, (b) Checkerboard (or 
plaquette) order 26,27 is related to stripes. While both 
can occur as stable phases (with similar energetics) in 
variants of the t-J model, 28,29 there are key experimental 
facts pointing towards stripe instead checkerboard order 
being the primary instability: (i) The STM data 9 clearly 
show that the rotation symmetry is locally broken from 
C4 down to C2. (h) The momentum-space pattern of 
spin excitations arising from checkerboard order has been 
shown to be incompatible with the neutron response of 
materials like La2-xBa x Cu04 or YBa 2 Cu306+,s- 20,30 



II. SUPERCONDUCTING VALENCE-BOND 
STATES 

Numerous microscopic calculations, for Hubbard or t-J 
models at low doping, have indicated a tendency towards 
states with broken translational symmetry. 18,31-36 At low 
temperatures T, it is conceivable that this coexists with 
superconductivity. 28 ' 36,37 In fact, such a scenario can be 
expected on general grounds: Upon destroying magnetic 
order in a square-lattice antifcrromagnet (AF), param- 
agnetic states with valence-bond (or spin-Peierls) order 
are known to appear. 17 The introduction of charge carri- 
ers by doping then leads to superconductivity, coexisting 
with bond order for a finite doping range. 36 ' 38 A global 
phase diagram has been worked out using a Sp(2A) 
mean-field theory applied to the t-J model supplemented 
by longer-range Coulomb interaction: 36 At small dop- 
ing, superconducting bond-centered stripe states occur, 
while homogeneous (Z-wave superconductivity is realized 
at larger doping. (Depending on microscopic parameters, 
the stripes may get replaced by checkerboard or spin- 
Peierls states.) Related superconducting charge-ordered 
states also appear in other theoretical treatments. 27,37 

Let us specify the various types of translational sym- 
metry breaking on a square lattice, assuming the mag- 
netic SU(2) symmetry to be unbroken, (a) Spin-Peierls 
states have a 2 x 1 unit cell where all sites are equivalent, 
but the bonds are modulated. The C4 rotation symmetry 
is broken, the ordering wavevector is Q = (n, 0),(0, it). 
(b) Stripe states have a N x 1 unit cell (N = 4 is par- 
ticularly stable), both sites and bonds arc modulated, 
Q = (±27r/A,0),(0,±27r/A), and C 4 is broken, (c) 
Checkerboard states have a N x N unit cell, both sites 
and bonds are modulated, Q = (±2n / N 7 ±2n / N) , but 
C4 is intact. While stripe and checkerboard states can in 
principle be either site- or bond-centered, experimental 
evidence points towards bond-centered structures. 9 ' 21,39 



In a quasiparticle picture, the symmetry-breaking 
orders can be translated into expectation values of 
fcrmionic bilincars: fa(\s) = (ck|C_kx) captures homoge- 
neous superconductivity, with fa (k) oc cos k x — cos k y for 
(Z-wave pairing, while fa(k) = (c k+ q a c\s.a) an d <fe(k) = 
(ck+Q,tc_ki) originate from charge order. The direc- 
tional dependence of fa (k) can be decomposed according 
to the representations of the point group. In the pres- 
ence of unidirectional order, i.e., C4 broken down to C 2 , 
the s-wave and (Z-wave order parameters inevitably mix. 
Our superconducting stripe states below will have 1/11,2,3 
all non-zero. (Note that e.g. fa alone generates a Fuldc- 
Ferrell-Larkin-Ovchinnikov state.) 



III. d-WAVE STRIPES 

Stripe states are best discussed in real space (phrased 
in the following for a square lattice of Cu atoms, keeping 
in mind that bonds of this lattice correspond to oxygen 
orbitals). Naively, the primary phenomenon of stripe 
order is a modulation of the on-site charge densities, 
(cJ CT Ci CT ), which translates into a momentum-independent 
(i.e. s-wave) order parameter fa (k). However, general 
arguments indicate 17,36 that the physics behind local or- 
dering acts on bonds instead of sites: Stripe formation 
is driven by the competition between kinetic and mag- 
netic energies, both living on lattice bonds. In such a 
bond-dominated stripe state, modulations in quantities 
like (c ia .Ci + A,cj} can locally have different signs on hori- 
zontal and vertical bonds, implying a (Z-wave component 
offa(k). 

Given the mixing of s- and cZ-wave components, stripe 
states may have ordering primarily in cither the s- or 
the (Z-wave channel. In the following, we shall argue 
in favor of stripes dominated by the (Z-wave component. 
Fourier-transforming fa(k) into real space 24 leads to a 
stripe state as in Fig. 1 - this is one of the main re- 
sults of this paper. Importantly, the ordering pattern in 
Fig. 1 appears perfectly consistent with the STM data 
of Kohsaka et al. (Fig. 4 of Ref. 9), where locally well- 
formed period-4 structures are seen, consisting of distinct 
ladder-like objects. 

Arguments in favor of stripes with a dominant (Z-wave 
component (Fig. 1) are: 40 (i) Energetics: In the hole- 
poor regions, strong horizontal bonds form due to the 
tendency towards dimerization (i.e. optimizing the mag- 
netic energy), whereas in the hole-rich regions strong ver- 
tical bonds form optimizing the hole kinetic energy, (ii) 
Coulomb interaction disfavors on-site charge modulations 
and thus suppresses the s-wave component. 

If stripe order coexists with superconductivity, then 
the bond pairing amplitudes, (cj-|-Cj+A,j.) 5 will be non- 
zero as well. Starting from a homogeneous (Z-wave super- 
conductor, stripe order will induce a modulated s-wave 
pairing component, described by ^(k). 
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FIG. 2: Doping dependence of quasiparticle gap energies from 
Sp(2iV) mean-field theory. The solid line shows the antinodal 
gap, i.e., the minimum QP energy along (k x ,iv), for t/J = 4, 
t' /t — —0.25 - here only spin-Peierls order, but no stripes 
occur for doping 5 > 3%. (For stripe states, the characteristic 
gap scale is more difficult to extract due to band folding.) In 
addition, the dashed line shows the minimum QP energy along 
the boundary of the AF Brillouin zone (BZ). 41 



IV. MEAN-FIELD THEORY 
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FIG. 3: Electronic structure of a 4 x 1 stripe state from mean- 
field theory, here with t/J — 1.95, t'/t = —0.25, overall dop- 
ing 8 = 1/8, and a hole distribution of 1.3 : 0.7 enforced by 
hand, see text, a) Dispersion along high-symmetry lines in 
k space, b) Constant-energy cuts through k space, showing 
nodal and antinodal structures. The spectra of horizontal and 
vertical stripes have been added, and an artificial broadening 
of J/20 was used. 



Superconducting stripe states can be obtained in 
Sp(2A r ) mean-field theory. 36 Consider an extended t-J 
Hamiltonian for fermions without double occupancies, 
Ci a , with spin a = 1 . . . 2N (N — 1 is the physical value): 
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with ri j = c\ a Ci a . The spin operators Si are fermion 
bilinears times the traceless generators of Sp(2iV). The 
fermion hopping, tij, will be non-zero for nearest and 
next-nearest neighbors, with values t and t' , while the 
exchange, Jy, is restricted to nearest- neighbor terms, J. 
The average doping 5 is fixed by J2ii n i) = NN S (1 — 5), 
where N s is the number of lattice sites. 

The spins and holes can be represented by auxiliary 
fermions fi a and spinless bosons bi, respectively, such 
that the physical electrons Cj Q = b\fi a . Via a Hubbard- 
Stratonovich decoupling of the AF interaction we intro- 
duce link fields Qij , defined on the bonds of the square 
lattice. After taking the limit N — ► oo, the slave bosons 
bi condense, (bi) — \/~Nbi, and the Qij take static saddle- 
point values. We are left with a bilinear Hamiltonian 
which can be diagonalized by a Bogoliubov transforma- 
tion. At the saddle point, the slave-boson amplitudes 
fulfill Y^i^f = N s 5, and the link fields are given by 
NQ tj = (J al3 fi a f]p), where J af3 is the antisymmet- 
ric Sp(2AT) tensor. The distribution of the bi is found 
be minimizing the saddle-point free energy, here the 
Coulomb repulsion Vij enters (on a classical level only). 

At non-zero doping and T = 0, all mean-field phases 
are superconducting. At larger t/J, the on-site charge 
distribution is homogeneous except for tiny doping, but 



bond-order of spin-Pcicrls type occurs at low doping. The 
doping dependence of the gap energy scale is shown in 
Fig. 2; it roughly follows the experimentally established 
pseudogap scale T*. 

At smaller t/J, stripe order occurs over a sizable dop- 
ing range. As in most mean-field theories, the tendency 
toward ordering is overestimated: In a large parame- 
ter regime, the mean-field theory predicts bond-centered 
stripes with maximal charge inhomogeneity and strong 
bond order. 36 To obtain the electronic structure for a 
more realistic stripe state from the mean-field theory, we 
have employed the following modifications: (i) The dis- 
tribution of on-site hole densities within the unit cell is 
enforced by hand through the bi values, (ii) The param- 
eter t/J is chosen such that the system is in the bond- 
ordered regime, but close to the transition to homoge- 
neous superconductivity. These modifications anticipate 
that both quantum effects and Coulomb repulsion (be- 
yond the classical approximation) reduce the amplitude 
of the charge modulations. 42 [Density-matrix renormal- 
ization group calculations for the 2d t-J model (with real- 
istic t') showed signatures for bond-centered stripes with 
reduced on-site modulation amplitude. 34 ' 35 ] 

A sample result for the electronic spectrum of a period- 
4 stripe state is shown in Fig. 3. In this calculation, the 
hole distribution was fairly inhomogeneous, b\ — b\ = 
1.3(5, b\ = b\ = 0.7(5, and the link fields \Qij\ varied by a 
factor of four within the unit cell. Thus, the stripe state 
contained both sizable s- and (i-wave components. Nev- 
ertheless, the spectrum is essentially gapless along the k- 
space diagonal, while a large gap appears near the antin- 
odes, Fig. 3. This toy calculation illustrates a key point, 
namely bond-centered stripe order is compatible with the 
presence of nodal quasiparticles, provided that the s-wave 
component of the order is not too large 43,44 - this com- 
patibility was subject of discussions in the past. 45 
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V. GLASSY DISORDER 

Let us comment on the influence of quenched disor- 
der arising from dopant impurities. In situations where a 
symmetry-breaking order is only "almost" static, adding 
disorder can pin the fluctuations and locally stabilize or- 
dered islands. The result is a state with static short- 
range order, not unlike in glassy systems. For stripe 
order, this means that domains with segments of hori- 
zontal and vertical stripes will coexist, with a checker- 
board structure of domain walls. 20,46 The electronic QP 
excitations now move in a static disorder potential, which 
couples strongly only to the QP near the antinodes, while 
nodal QP are little affected. As a result, the scattering 
rate along the normal-state Fermi surface will be strongly 
energy-dependent, and gapless coherent nodal QP will 
coexist with incoherent antinodal QP. 

Importantly, the disorder potential acts as a ran- 
dom field on charge stripes, smearing out any finite- 
temperature phase transition in the charge sector. 18 ' 46 
(This is distinct from the spin sector, where static order 
at low T is still accompanied by a sharp phase transi- 
tion, as the disorder is of random-mass type.) At ele- 
vated temperatures, signatures of local valence-bond (i.e. 
singlet) formation will be visible below a temperature 
T* ~ J, and will likely evolve continuously from the high- 
temperature "pseudogap" scale to the low-temperature 
"superconducting" gap. 

VI. THEORY VS. EXPERIMENT 

The described picture of valence-bond order, coex- 
isting with superconductivity at low T, nicely ties in 
with various features of recent experimental data, (i) 
Models of coupled spin ladders, arising from bond- 
centered stripes, 47,48 provide an excellent description of 
the spin dynamics in stripe-ordered Lai 5 / 8 Ba 1 / 8 Cu04. 21 
(ii) The short-range charge order seen in STM on un- 



derdoped P^S^CaC^Os+a and Ca2- 2: Na a ;Cu02Cl2 9 is 
bond-centered and locally breaks the lattice rotation 
symmetry down to C2 The glassy real-space structure 
is compatible with impurity-pinned short-range order, 
where fluctuations in the charge order parameter are pri- 
marily of phase instead of amplitude type. 20 (iii) The 
apparent d-wave gap of the non-superconducting stripe 
compound La 15 / 8 Ba 1 / 8 Cu04 is naturally explained by 
d-wave stripes, (iv) To explain the absence of 3d super- 
conductivity in La 15 / 8 Bai/8Cu04, "antiphase supercon- 
ductivity" has been proposed. 49 This is similar to 03 (k) 
above, and broadly consistent with the picture advocated 
here. 



VII. CONCLUSIONS 

Building on earlier work on charge order in cuprates, 
we have pointed out that bond-ordered stripe states, pos- 
sibly of glassy character and coexisting with supercon- 
ductivity, provide a phenomenological framework which 
is consistent with many features of recent cuprate ex- 
periments. We believe that this strengthens the case for 
bond order (and associated local singlet formation) to 
be a common tendency in cuprates, likely relevant also 
for pseudogap phenomena. Let us note, however, that 
we think that stripe-like translational symmetry break- 
ing is a phenomenon competing with superconductivity, 
i.e., pairing is suppressed by stripes with spatial long- 
range order. 
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